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Abstract. - We propose a straightforward and efficient procedure to perform dynamical 
mean-field (DMFT) calculations on the top of the static mean-field LDA+U approximation. 
Starting from self-consistent LDA+U ground state we included multiplet transitions using the 
Hubbard-I approximation, which yields a very good agreement with experimental photoelectron 
spectra of <5-Pu, Am, and their selected compounds. 

Introduction. - Electronic, magnetic and superconducting properties of actinide elements 
recently attracted significant interest and attention in the condensed matter physics. Most 
intriguing are the phenomena at the localization threshold of the 5/ series, which is crossed 
between Pu and Am, where the electron-electron correlations play a prominent role [1,2]. 

During last few years, electronic structure calculations of Pu and Am based on the conven- 
tional band theoretical methods (the local density or generalized gradient expantion approxi- 
mations LDA/GGA to the density functional theory) could not explain essental experimental 
data. While, the LDA/GGA band structure calculations predict a local magnetic moment 
(ordered or disordered) to form on the Pu [3] and Am [4] atoms, none of them were seen in 
the experiment [5]. Also, the same papers attempted to evaluate the photoemission spectra 
(PES) and electronic specific heat in Pu and Am making use of single-particle LDA/GGA den- 
sities of states (DOS), incorrectly assuming weak electron correlation character of 5/ systems 
at the borderline between the localized, nonbonding, behaviour and the bonding situation of 
electronic bands. 

It was shown recently that the around-mean-field (AMF)-LSDA+U correlated band theory 
gives non- magnetic ground state for Pu [6] and Am [7] . Also, the equilibrium volumes and bulk 
moduli for S-Pvl and /cc-Am are calculated in a good agreement with experiment. However, 
there is a clear disagreement between the AMF-LDA+U calculated DOS and PES, questioning 
the validity of this approximation. The DMFT calculations [8] using the (AMF)-LSDA+U 
form of interacting Hamiltonian demonstrated that experimental PES and high 7 coefficient 
of the electronic specific heat in S-Pu and its selected compounds [9, 10] originate from the 
excitations, and not from the ground state DOS alone. While accounting for dynamical 
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fluctuations, the DMFT of Ref. [8] does not take into account atomic-like excitations which 
can play an important role in Pu and Am. 

In this Letter, we develop a DMFT based computational scheme based on multi-orbital 
Hubbard-I approximation (HIA) [11-13] including the spin-orbit coupling (SOC) which ex- 
plicitly accounts for the atomic-like multiplet transition excitations in Pu, Am and their com- 
pounds. Starting from the non-magnetic ground state calculated with the static mean-field 
AMF-LSDA+U approximation, we obtain excitation spectra of Pu and Am in surprisingly 
good agreement with PES, in support of the atomic-like origin of the electronic excitations in 
these materials. 

Methodology. - We start with the multi-band Hubbard Hamiltonian [11] H = H° + H lnt , 
where 

H ° = E E H kiP c UiP = E E ff°,,(k)4(kMk) , (1) 

i,j a.0 k a.0 

is the one-particle Hamiltonian found from ab initio electronic structure calculations for a 
periodic crystal. The indices i,j label the lattice sites, a = (imcr) denote the spinorbitals, 
and k is the k-vector from the first Brillouin zone. It is assumed that the electron-electron 
correlations between s, p, and d electrons are well described within the density functional 
theory, while the correlations between the f electrons have to be considered separately by 
introducing the interaction Hamiltonian 

^ mt = ^E E ^( m ^ m 2\Vr\m 3 ,m 4 )cl mia cl m2(yl c lm4lj/ c lm3a . (2) 

i mi ,m2, 7713,7714 (7,<7' 

The V ee is an effective on-site Coulomb interaction [11] expressed in terms of the Slater 
integrals Fk and the spherical harmonic \lm). The corresponding one-particle Green function 

G(k,z) = (z + //-if°(k)-E(k,z)) 1 (3) 

is expressed via H° and the one-particle selfenergy S(k, z) which contains the electron-electron 
correlations, where z is a (complex) energy with respect to the chemical potential /i. The 
interactions J2J act only in the subspace of f-states. Consequently, the selfenergy S(k, z) is 
nonzero only in the subspace of the f-states. 

The simplest mean-field approximation (often called L(S)DA+U) [14] neglects the k- and 
energy-dependence of £ replacing it by the on-site potential V+u- For a given set of spin- 
orbitals \ma) the potential reads: 

[V+uY mm > = E ((™,P\V ee \m',q) ~ (m,p\V ee \q,m')S^)n;[ q , (4) 

p-q-o' 

where n„ m / is the local orbital occupation matrix of the orbitals \mcr). It was shown in 
Ref. [15] that the Kohn-Sham equation with the potential Eq.(@J can be obtained by making 
use of variational minimization of the LDA+U total energy functional (i.e. of the expectation 
value of the multiband Hubbard Hamiltonian) in a way similar to the conventional density 
functional theory [16]. 

In what follows we use the local approximation for the selfenergy, i.e., we assume that it is 
site-diagonal and therefore independent of k. Then we can employ the "impurity" method of 
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Ref. [11]. We first obtain a local Green function integrating G(k, z), Eq.Q, over the Brillouin 
zone 

G(z) = -^— f dk(W M -#°(k) ~* (5) 
Vbz Jbz v ' 

and define a "bath" Green function (the so-called Weiss field) Qo(z) 

g (z)=(G- 1 (z) + E(zj)~ 1 . (6) 

All Q/, G, and S in Eqs. (|S1EJ) are matrices in the subspace of the \ma) f-orbitals. 

The DMFT self-consistency condition is now formulated by equating G(z), Eq.lJSJ, to the 
Green function G(z) of a single-impurity Anderson model (SIAM) [17] 

G(z) = (z + ^-eo- A(z)-S(z)) _1 (7) 

where A(z) is the effective hybridization function and is the SIAM self-energy. We write 
the Eq.0 in the form of the Eq.© 

(z + M-eo-A(z)) _1 =Go(z)= (g- 1 (z) + S(z)) _1 (8) 

from which it follows that Qo(z) has the meaning of the SIAM "bath" Green function. The 
iterative procedure to solve the periodic lattice problem in the DMFT approximation is now 
formulated in a usual way [11]: starting with single particle Hamiltonian H° (k) and a guess 
for local £, the local Green function is calculated using the Eq.|0 and the "bath" Green 
function is calculated from Eq.©; the SIAM is solved for this "bath" and a new local £ is 
calculated from Eq.lJTJ), which is inserted back into Eq.j5j. 

The LDA+U procedure can be viewed in the same way. There is no need to apply the 
full DMFT iterative procedure described above and to solve the SIAM with the "bath" from 
Eq.JHJ). The self-energy T,(z) is now approximated by a static potential V+u from Eq.Q, 
and the well-known relation between the Green function, Eq.@), and the local orbital oc- 
cupation matrix n a mm , = — 7r _1 Im J M dE G(E)° lm , is used. In addition, the charge- and 
spin-densities needed to construct the single-particle Hamiltonian Ho(k) in Eq.jsj are calcu- 
lated self-consistently. We emphasize that LDA+U approximation is generically connected 
with the LDA+DMFT procedure. 

Hubbard-I approximation. - Here we attempt to build a computational scheme based on 
self-consistent static mean-field LDA+U ground state that will allow us to access the correlated 
electron excitations. We specifically choose the around-mean-field version of LSDA+U (AMF- 
LDA+U) which was shown to describe correctly the non-magnetic ground state properties of 
<5-Pu [6], fee- Am, Pu-Am alloys [7], and selected Pu-compounds [9]. We extend our previous 
works [6,7] towards the DMFT to account for the multiplet transitions which are necessary for 
a correct description of PE excitation spectra. We use the multiorbital HIA which is suitable 
for incorporating the multiplet transitions in the electronic structure, as it is explicitly based 
on the exact diagonalization of an isolated atomic-like f-shell. 

Further, we restrict our formulation to the paramagnetic phase, and we closely follow the 
procedure described in [11]. We construct the atomic Hamiltonian including the spin-orbit 
coupling (SOC): 

er,<r a,a 
i jat = £(l-s)m 1 C C m 1 ff C m2 «r' + ( m l m 2 l m 3™4) c L 1( xcj„ 20 ., C mia , C m3(T , (9) 

mi -,7*12 mi. . .7714 
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and perform exact diagonalization of H &t \v) = E„\v) to obtain all possible eigenvalues E v 
and eigenvectors \v). 

The HIA "chemical potential" fin is then calculated as 

(n) = -Tr [N exp(-/?[ff at - /«-#])] (10) 

for a given number of particles (n). Here, (3 is the inverse temperature and Z is the partition 
function. 

Finally, the atomic Green function is calculated as follows: 

* ' 1 («|c m , a \v) McL „i\u) 

l 1 ^ \ z )\m im 2\ Z ) 



Z Z+(E^- /iff A^) - [E v - fi H N v ) 



x [exp(-/3(E„ - hhNu)) + exp(-/?(£ M - JV M ))] (11) 
and the atomic self-energy is evaluated as: 

[£*(*)] W = ^m ima *^ ~ C(l • s)™^ - [(G at (z))"T " . (12) 

L \ /J mim2 

Assuming that the self-consistent LDA+U calculations are performed, the LDA+U Green 
function is evaluated as 

G+u(z) = t^[ dkfz + fi-HoM-V+uY 1 . (13) 

VBZ JBZ K ' 

In Ea. l|13|) we took into account the presence of SOC for both Hq(\s.) and V+jj as described 
in Ref. [6]. We used LDA+U eigenvalues and eigenfunctions calculated in the full-potential 
LAPW basis [6,15] to construct the on-site spin-orbital Green function matrix, Ea. (|13fl . 

We evaluate the static "bath" Go(z) from Eq.® and find the hybridization function A(z) 
together with (eo — /i) which determines the energy of "impurity" level with respect to the 
solid potential: 

(z + ii - e +u - e(l • s) - A(z)) ^ - g%{z) = (G-]j{z) + V +U (zj) ^ . (14) 

Here we added the SOC explicitly and assumed the paramagnetic case. 

We point on a difference in a physical meaning of eo m the original SIAM and in the 
auxiliary SIAM used in LDA+DMFT (LDA+U): in the former, it labels the position of the 
non-interacting f(d)-level so that the chemical potential is determined self-consistently for a 
given A(z); in the latter, the chemical potential fi is determined by a periodical crystal - 
basically by the Green function of Eq.© - so that the auxiliary eo accommodates all the 
electrostatic shifts between correlated electrons and the potential of the solid. 

Now we can formulate a simple approximate procedure to solve the DMFT Eas. (|5l6l7|l 
with Hubbard-I self-energy, Eq. l|12|) . We assume that the self-consistent static mean-field 
LDA+U already gives a correct number of particles and hybridization. Using the LDA+U 
Green function, Eq.©, and the potential, Eq.(@J, we evaluate the Weiss field Q^_ u (z). Then 
we insert the HIA self-energy calculated for the same number of correlated electrons as given 
by LDA+U, into this "bath" , and calculate the new Green function 

G(z) = (iG^iz)}- 1 + (e +u - e H ) E H (z)) ' , (15) 

where (e+u — ejy) is chosen so as to ensure that n = 7r _1 Im J 11 dETr[G(E)} is equal to a 
given number of correlated electrons [19]. 
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Valence band photoemission spectra. Starting from self-consistent AMF-LSDA+U 
ground state solutions for 5-Pu [6], and fee- Am [7], and making use of the corresponding 
eigenvalues and eigenfunctions we evaluate the G+u(z) Green function given by Eq. l|13|l and 
the Go U {z) "bath" Green Function (Ea.fity) [20]. In HIA calculations Eq.0- EJ we used 
the commonly accepted values of SOC constants £ = 0.3 eV (Pu) and 0.34 eV (Am). To find 
Hh Ea. flOp . we choose the AMF-LDA+U values of (n), namely n/=6.0 for Am and n/=5.4 
for Pu [21]. The self-energy Ea. (|12fl was calculated along the real axis for z = E — E-p + iS, 
where 5 = 63 meV [18], and (3 was varied from 100 to 1000 eV -1 . We found no sizable effect 
due to the variation of f3 for a given in the resulting spectral density Ea. i|15fl . 

In Fig. la we show the f-projected DOS (fDOS) from AMF-LDA+U together with spectral 
density calculated from Ea. H15|) . For fee- Am, the well localised fDOS peak at -4 eV transforms 
to the multiplet of excited state transitions / 6 — > / 5 below the Fermi energy and fDOS- 
manifold around +2 eV to / 6 — > f 7 multiplet transitions. Although there is no doubt that 
the 5f multiplets must dominate the experimental valence-band spectra of Am-based systems, 
individual lines are not resolved (except for partly resolved features in the spectrum of Am 
metal), and the position of the 5f intensity in the energy spectrum is the main indicator of 
the agreement with calculations. In this sense, the calculated Am spectral density is in a 
good agreement with PES [22]. It also agrees with similar calculations [18] [23]. Furthermore, 
we performed the calculations of AmN and AmSb using the same {U, J} set of values as for 
elemental Am. The AMF-LDA+U yields AmN as an indirect gap semiconductor and AmSb 
as a semi-metal. The HIA spectral densities for Am / manifolds in AmN and AmSb are shown 
in Fig. lb, c and are in good agreement with PES experiments [22]. 

Experimental valence-band spectra of <5-Pu and several other Pu systems exhibit three 
narrow features within 1 eV below Ep, the most distinct one very close to Ep being accom- 
panied by a weaker feature at 0.5 eV and another one at 0.8-0.9 eV. Their general occurrence 
and invariability of characteristic energies practically excludes any relation to individual fea- 
tures in density of electronic states. Instead, a relation to final state multiplets has been 
suggested [24,25] among other possible explanations . Similar to Am, the link to atomic mul- 
tiplet is corroborated by the present calculations also for Pu. For <5-Pu shown in Fig. 2a, the 
AMF-LDA+U fDOS manifold at around -1 eV transforms into a set of multiplet transitions 
with high value of spectral density at Ep. Similarly to the Pu metal, the PES exhibits the 
three most intense peaks for a broad class of Pu compounds. The calculations performed for 
PuTe (see Fig. 2b) indeed demonstrate the three-peak pattern similar to Pu, in agreement 
with experiment [26] . Also we point out a good agreement between our calculations and recent 
DMFT calculations of Ref. [8], as well as those of Svane for PuSe [18]. One should note that 
the high spectral intensity at the Fermi level is explaining enhanced values of the 7-coefficient 
of electronic specific heat, observed in S-Pu and other Pu systems. 

As a conclusion, we have shown that the three narrow features observed in the valence band 
spectra of Pu and majority of Pu compounds can be identified with the most intense atomic 
excitations (multiplets), calculated using the LDA+U and the multi-orbital HIA calculations. 
The calculations explain that the atomic excitations can be observed even if the 5f states are 
not fully localized as in <5-Pu, and the atomic character fixes the characteristic energies (not 
intensities) such that similar features are found in spectra of diverse Pu systems. A reasonable 
agreement with experiment is found also for Am and its compounds calculated on the same 
footing. 
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Fig. 1 - Color-on-line. The fDOS (AMF-LDA+U) and Spectral Density (HIA) for fee- Am (top-a), 
AmN (middle-b) , AmSb(bottom-c) calculated along z = E — Ef + i8 where 5 = 63 meV is used. 
Experimental PES (arb. units) from Ref. [22] are shown [27]. 

Fig. 2 - Color-on-line. The fDOS (AMF-LDA+U) and Spectral Density (HIA) for 8-Pu (top-a) and 



